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lO ■ The low number density of the Sloan Digital Sky Survey (SDSS) Lumi- 

^ ■ nous Red Galaxies (LRGs) suggests that LRGs occupying the same dark mat- 

■ ter halo can be separated from pairs occupying distinct dark matter halos with 
high fidelity. We present a new technique, Counts-in-Cylinders (CiC), to con- 
strain the parameters of the satellite contribution to the LRG Halo-Occupation 

■ Distribution (HOD). For a fiber collision-corrected SDSS spectroscopic LRG 
j_] ■ subsample at 0.16 < z < 0.36, we find the CiC multiplicity function is fit 

■ — ' by a halo model where the average number of satellites in a halo of mass 

M is {N,at{M)) = ((M - Me„0/Mi)" with M,^t = 5.0+JJ(+|^) x IO'^Mq, 
Ml = 4:.95tofe{toil) x W^Mq, and a = 1.035l°:J°(l°:^t) at the 68% and 95% 
confidence levels using a WMAP3 cosmology and z = 0.2 halo catalog. 

Our method tightly constrains the fraction of LRGs that are satellite galax- 
ies, 6.36l°:|%, and the combination M^^t/IO^^M© + a = 1.53+°;°^ at the 95% 
confidence level. We also find that mocks based on a halo catalog produced 
by a spherical overdensity (SO) finder reproduce both the measured CiC multi- 
plicity function and the projected correlation function, while mocks based on a 
Friends-of- Friends (FoF) halo catalog has a deficit of close pairs at ~ 1 Mpc/h 
separations. Because the CiC method relies on higher order statistics of close 
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pairs, it is robust to the choice of halo finder. In a companion paper we will 
apply this technique to optimize Finger-of-God (FOG) compression to eliminate 
the 1-halo contribution to the LRG power spectrum. 



Introduction 



The Sloan Digital Sky Survey (SDSS; I York et al.l (120001 )) has recorded the largest sam- 



pie of Luminous Red G alaxies (LRGs), probing a volume of ~ 1 (/;, ^ Gpc)'^ out to 2; ~ 0.5 



(lEisenstein et al.ll200ll ) and making it ideal for studying large scale structure. Understanding 
the small-scale relationship between the galaxy and dark matter density fields is essential to 
extracting t he linear matter power spectrum from the ga laxy power spectrum, even on very 
large scales JSchulz fc WhitJbood : ISanchez fc Col^boosh . 



The Hal o Occupation Distribution (HOD) is a popular and useful description of this 
relationship (ISeljakl I2OOOI : iPeacock fc SmithI I2OOOI : ICooray fc ShethI |2002| ) , and can be used 
to constrain the rat e of merging, d i srupt ion or evolution in well-d efined galaxy populations. 
Conrov et all ((2003), Iwhite et aP ((2003), and IWake et al.l ((20081) have used this framework 
to constrain the rate at which LRGs merge or are disrupted in clusters. iBrown et al.l (120081 ) 
combine HOD constraints with luminosity function measurements of red galaxies to deduce 
that stellar mass build-up in clusters occurs primarily in the satellite galaxies or intraclus- 
ter light betw een z = 1 and z = 0, while the central galaxies grow only modestly, with 



L,en ~ Ml/3. 

and SDSS 



Zheng et al.l (120071 ) constrain stellar mass gro wth between pEEP2 



[z ~ 0) galaxies using the HOD description, and IConroy et al.l (120081 ) employ 
HO D mod e hng t o illu minate the fate of 2 ~ 2 star-forming galaxies. Others researchers. 



ChenI (l2007h and IHo et a. 



(120071). use the HOD description to study the spatial dis- 



Wake et al.l (120081 ) also argue that the small-scale clustering 



e.g. 

tribution of satellite galaxies, 
at different redshifts constrains the scatter in halo merger histories which can be compared 
with predictions of hierarchical models. 



S everal groups have used two and three point statistics (iBlake et al.ll2008l: iKulkarni et al. 



20071 : IWhite et al.ll2007l : IZheng et al 



well as galaxy- galaxy lensing (IMandelbaum et al. 



2008 



Wake et al. 



2008 



Padmanabhan et al 



200 



as 



20061 ) to constrain the HOD of LRGs. 



Ho et al.l (120071 ) have taken a more direct approach and used X-ray determined cluster masses 
to measure Nlrg{^)- Though these analyses were performed on samples with different lu- 
minosity and redshift ranges, they offer seemingly conflicting results on both the slope a at 
the high mass limit of the satellite term Nsat ~ M° and the fraction of LRGs that are satellite 
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galaxies. IHo et al.l (120071) find a 



Kulkarni et al 



mm find a 



mm find ^ 1.8 f or a« 



0.6 when fit t ing th e total LRG number Nt miM) oc M' 
J2008h find a ~ 2.1 - 2.6, and lzheng eVal. 



1.4, 


Blake et al. 


Kulkarni et al. 



(120071 ) rep ort a satellite fracti on of ~ 17% 



while iBlake et all (l2008f )'s redshift slices span 3-8%, and IZheng et al.l mm find 5 - 6% 



for the LRG subsample studied in this paper. The most luminous elliptical galaxies in 



Mandelbaum et al.l (120061 ) have a satellite fraction of < 10%. 



The low number density of the SDSS Luminous Red Galaxy (LRG) sample suggests 
that LRG pairs occupying the same dark matter halo can be separated from pairs occupying 
distinct dark matter halos with high fidelity. In this paper we explore that intuition, and 
show that one-halo pairs can be identified with ~ 75% completeness and < 27% contamina- 
tion by simple cuts in the transverse separation Ar_L and LOS separation Ary. Furthermore, 
these pairs can be grouped together using a Friends-of-Friends (FoF) algorithm to estimate 
the LRG group multiplicity function. We apply this technique to a sample of LRGs from 
SDSS to constrain their HOD. We find that both the high values of a ~ 2 and high satellite 
fractions reported in previous papers are inconsistent with the 0.16 < z < 0.36 SDSS LRG 
group multiplicity function measured here. In contras t to previous rn e thods which rely on 
2 and 3 point statistics to constrain the HOD (as in iKulkarni et al.l (120071 )). our method 
probes the HOD more directly by estimating the group multiplicity function from the higher 
order statistics in the LRG density field in the one-halo dominant regime. 



We present an overview of the CiC method in § 12.11 and apply it to an approximately 
volume limited subsample of SDSS LRGs in § 12.21 addressing the comphcations of fiber 
collisions, incompleteness, and complex angular masks. The CiC technique developed here 
requires calibration on mock galaxy catalogs. We summarize our A^-body simulation pa- 
rameters in § 12. 3[ § 12.41 presents the HOD model we employ throughout this analysis and 
details how we populate our simulations with galaxies. § 12.51 describes the CiC technique 
to measure the LRG group multiplicity function and its calibration with simulations. The 
HOD parameters are fit using a maximum likelihood analysis explicated in § 12.61 In § [3] we 
present the CiC multiplicity function of our SDSS LRG subsample and describe the relation 
between the CiC and true group multiplicity functions. We present the constraints on the 
HOD parameters and their implications for the fraction of LRGs that are satellites, as well 
as the mass distribution of halos hosting LRG groups with Usat satellites. Mock catalogs 
produced using the C iC maximum l i keliho od HOD and a spherical overdensity (SO) halo 
catalog agree with the lMasjedi et al.l (120061 ) measurement of Wp{rp) for this sample when the 
large scale bias is adjusted with a single parameter for the central galaxy HOD. In § 13.51 we 
compare these results with a mock LRG catalog based on a FoF halo catalog and with other 
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HOD measurements in the literature. We show that while the FoF and SO catalogs can 
both match the observed CiC multiplicity function, the FoF catalog produces mock catalogs 
with a deficit of halos at 1 Mpc//i that is evident in the projected correlation function. In 
§ m we comment on the strengths and weaknesses of the CiC method and summarize our 
conclusions in § [5l 



Throughout this paper we adopt the ISpergel et al.l (120071 ) cosmological parameters used 
in our simulations to convert redshifts to distances: {flrn,^b,^A,ns, as, h) = (0.26, 0.044, 
0.74, 0.95, 0.77, 0.72). All distances and separations are in comoving coordinates. 



2. Methodology 

2.1. Overview of the Method 

The goal of this section is to measure the group multiplicity function of a subsample of 
the SDSS LRGs. For a complete spectroscopic sample covering the full sky (or in a periodic 
simulation box), our method is as follows: 



• Identify 'one-halo' pairs of galaxies. In what follows, our criteria for two galaxies to be 
a one-halo pair is Ar_L < 0.8 Mpc/h and Ary = 20 Mpc/h (equivalently Az/(1 + z) < 
Pmax = 0.006), where both are comoving separations. These choices were motivated 
by results on mock LRG catalogs and will be discussed in more detail in later sections. 

• Group pairs of galaxies into groups using a FoF algorithm. The number of groups with 
Usat satellites, Ncic{nsat) is the CiC multiplicity function. 



In § 12. 21 we present the technical details of accounting for the facts that the SDSS has bound- 
aries and holes, that the spectroscopic sample of LRGs is incomplete, and that the SDSS 
cannot simultaneously take spectra of two objects separated by < 55", so that regions of the 
sky observed only once spectroscopically may have missing close pairs of LR Gs. We use the 



LRGs from the SDSS imaging sample to supplement the spectroscopic sample (IBlanton et al. 



20051 : lAdelman-McCarthy et al.l 120071 ). We identify potential pairs from the imaging sample 



and calibrate this step using pairs of objects from the spectroscopic sample. Since most non- 
isolated LRGs are in groups of 2 and candidate pairs from the imaging sample neighboring 
more than one LRG are highly likely to be group members, we apply the small correction for 
false LRG pair detections to Ncic{f^sat = !)• A more complex scheme involving corrections 
at each Usat would not have enough statistics to calibrate on the spectroscopic sample. We 
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eliminate from our sample LRGs close to the survey boundary, though they are allowed to 
be grouped with LRGs away from the boundary. This ensures that our multiplicity function 
is not biased due to unobserved LRGs outside the boundary. However, since the bright star 
masks are numerous and individually very small, this approach is not practical for dealing 
with objects near bright star masks. Instead, we adjust N{nsat) by estimating the probabil- 
ity that there is an LRG covered by each bright star mask, and then computing the change 
in N{nsat) if there were one. Table [1] shows that all of these corrections are small. 



2.2. Data 



The SDSS (IStoughton et al.ll2002l : lAbazajian et al.ll2003l ) has imaged ~ 10^ deg^ in u 



g, r, i, and z. From th is sample, spectroscopi c LRG targets are efficiently selected using two 
color /ma gnitude cuts (jEisenstein et al.ll200ll ). The tiling algorithm ensures nearly complete 
samples ( Blanton et all 2003 ). However, spectroscopic fiber collisions prohibit simultaneous 
spectr oscopy for objects se parated by < 55", leaving ~ 7% of targeted objects without red- 
shifts ( iMasjedi et al.ll2006l ). Overlapping plates on ~ 1/3 of the survey area mitigate this 
problem and permit us to calibrate this effect in our analysis, as detailed below. The 'pho- 
tometric sample' as referred to below consists of objects from the i maging sample that were 
targeted as LRGs according to the color /magnitude cuts laid out in lEisenstein et al.l (120011 ) 
but lack spectra. The 'spectroscopic sample' consists of objects from the imaging sample 
that were targeted as LRGs and subsequently observed. 



The goal of this analysis is to measure the group multiplicity function for the spec- 
troscopic LRG sample with —23.2 < Mg < —21.2 and 0.16 < z < .36. This sanaple is 
app roximately vol u me-lim ited with n ^ 9.7 x 10^^[h^^ Mpc)^^, and IZehavi et al.l (120051 ) 
and iMasjedi et al.l (120061 ) have measured the projected correlation function for this sample 
on small and intermediate scales. We begin with the entire set of LRG t arget galaxies for the 



DR4+ sample from the NYU V alue-Added Galaxy Catalog (VAGC) (IBlanton et al.l 12005 



Adelman-McCarthy et al.l 120071 ) so that our sample includes the complete sample of galax- 



ies satisfying our Mg and z cuts. We correct our multiplicity function statistically for the 
inclusion of a small number of interlopers. 



Our analysis requires the identification of all close pairs of LRGs in the sample, where 
close pairs satisfy A9 < 6max{,z) and Ary < Ary max- The dominant source of LRG close pair 
incompleteness is from LRG-LRG fiber collisions; this results in a deficit of pairs separated 
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by < 55" in the spectroscopic sample. LRG-main galaxy fiber collisions or spectroscopic 
incompleteness result in further pair incompleteness. Finally, the bright star mask occupies 
1.88% of the survey area. We statistically correct for the unobserved LRGs in these regions. 



Following lEisenstein et al.l (120051 ) we cut sectors with < 60% spectroscopic complete- 
ness. These sectors lie primarily along the boundary of the survey, and remove only 276 
objects from our sample. The remaining survey area is 5564 deg^, of which 104 deg^ is 
covered by a bright star mask. Using inverse random catalogs from the VAGC to trace the 
survey geometry, we remove all objects within 9max{zmin = 0.16) of the survey boundary, 
which amounts to 469 deg^, or 8.6% of the total area. However, we keep track of all LRGs 
in the boundary, because they are allowed to form pairs with the objects remaining in our 
sample and therefore contribute to the group multiplicity function. Only 53 objects from 
the boundary were included in CiC groups. Note that while there are many corrections for 
incompleteness, they are all small and well-constrained. 



2.2.1. Sample Model and Color /Magnitude cuts 



There are 41721 objects in our spectroscopic sample passing the redshift, M^, sector 
completeness, and boundary cuts, and 8167 objects from the photometric sample passing 
the same sector completeness and boundary cuts. To reduce the level of contamination 
in our LRG multiplicity function, we apply color and magnitude cuts to objects from the 
photometric sample falling within A6max of a spectroscopic LRG to select a high-fidelity 
sample of close pair candidates. Sin ce we know the redsh ift of the proposed group from 
the spectroscopic neighbor, we use cy (lEisenstein et al.ll200ll ) as a redshift indicator and r^et 
of the photometric sample object as an absolute magnitude indicator at the spectroscopic 
redshift. 

C|| =0.7(^-r) + 1.2(r-i-0.18) (1) 

Using the colors of objects in our spectroscopic sample, we find the upper and lower limits of 
C|| that encompass 95% of the spectroscopic LRGs as function of redshift in bins of Az = 0.01. 
This relation is shown in Fig. [H For a candidate close pair match, we discard photometric 
sample objects with cy falling outside this region, gi ven the redshift o f the spectroscopic 
group member. Since the k-|-e corrections described in lEisenstein et al.l (l200ll ) have already 
been applied to the spectroscopic sample, we estimate the expected Mg of a photometric 
object at redshift z of the spectroscopic group member by 



M, 



g, photo 



■^g,spec ~l~ (j'petjphoto ^pet,spec) 



(2) 
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Fig. 1. — The mean cy — z relation and bands including 95% of objects in each = 0.01 
bin. C|| (Eqn. [1]) is used as a redshift indicator for objects targeted as LRGs but lacking 
spectra. 
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and disregard objects with Mg^photo < —23.2 or Mg^photo > —21.2. 



For a random sample of LRG targets, 36.6% of them will pass the lEisenstein et al. 



(120011 ) LRG cuts as well as our sample's magnitude and redshift cuts, based on the ratio 
in DR4+ for targeted LRGs that have spectra. The fraction of objects in our photometric 
sample that would pass our spectroscopic magnitude and redshift cuts will be slightly higher 
than for the random sample of LRG targets, since it includes fiber collision pairs very likely 
to be LRGs. We compute this fraction for the photometric sample at the end of the group- 
finding algorithm to circumvent this issue. We model the photometric sample as composed 
of Npass objects which would pass the absolute magnitude and redshift cuts of our sample, if 
they had spectra; the remaining objects we label as Nfaii- A large fraction of objects in the 
Nfaii group are below our minimum redshift cut Zmin = 0.16. Tests on targeted LRGs with 
spectra that fail our LRG subsample cuts show that the Nfau group is well-approximated as 
uncorrelated with the spectroscopic LRG sample. 



Candidate close pairs of spectroscopic LRGs from the photometric sample are naturally 
divided into two groups. The first, denoted 'FB' for 'fiber', includes objects 9 < 55" from 
a spectroscopic object. These objects do not have spectra primarily due to a fiber collision 
with the neighboring spectroscopic LRG. The bulk of candidate close pairs fall into this 
category, and the contamination for such pairs by objects at different redshifts is low. The 
remaining close pairs have separations 55" < 9 < 9max{zspec), arise from the small overall 
incompleteness of the survey or fiber collisions with MAIN galaxies, and are denoted 'INC. 
The probability of contamination is larger for this type, but still manageable. We use all 
targeted LRGs with spectra but failing our Mg or z cuts to compute an average rate at 
which an Nfau object will pass our color /magnitude cut on the photometric sample by 
comparison with the spectroscopic LRG pairs in each type of collision: Ppass,FB = 0.078 and 
PpassjNC = 0.070. The rate is slightly lower for INC collisions because they are more likely 
to be at a lower redshift where cy is more discriminating, since 9max{z) decreases with z. The 
number of close pair contaminants we expect from this sample for our two types of collisions 
are 

NpEJail = . * Nfail *Ppass ~ 2.1 (3) 

-^survey 

NiNCJail = — ^ * Nfail * Ppass ~ 26.8 (4) 

^survey 

where ApB and A^jig are the total areas enclosed by annuli of 55" and 9max{zspec) around 
each LRG, and Asurvey is the total survey area after removal of the boundary. Nfau was 
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estimated at the end of the group-finding algorithm; see below. 

To estimate the completeness and contamination for the Npass sample, we rely on ob- 
served collisions from the spectroscopic sample. For FB collisions, we found 295.5 colli- 
sions between objects in the spectroscopic sample that passed the color/magnitude cuts. 
276 of these passed the Az cut, so 93.4% of objects satisfying A6 < 55" and passing the 
color /magnitude cuts would be considered one-halo pairs if the redshift were known. Fur- 
thermore, 19 pairs failed the color /magnitude cuts but passed the Az cut. Therefore, the 
sample is 93.5% complete after the color /magnitude cuts and 6.6% contaminated. For INC 
collisions, we found 2494 collisions passing the photometric color/magnitude cuts; 2020 of 
these passed the Az cuts, so would be assigned as pairs. 110.5 pairs passing the Az cuts 
failed the color magnitude cuts. Therefore, this collision sample should be 94.8% complete 
and 19.0% contaminated. 

Our approach is to estimate the number of interlopers and apply the correction to the 
number of groups of 2 LRGs. 568 objects from the photometric sample are < 55" from a spec- 
troscopic LRG object and pass the color /magnitude cuts at the spectroscopic redshift; 250 
objects from the photometric sample are 55" < 9 < 9max{z) and pass the color /magnitude 
cuts at the spectroscopic redshift. Based on the rates measured from spectroscopic collisions 
and after correcting for the expected number of interlopers NpBjaii + Nijqcjaiu "we expect 
709.3 of the remaining 789.1 collisions to be 'true' collisions (i.e., would pass the Az criterion 
as well, if the redshift were measured), and we expect to have missed 46.3 true collisions due 
to our color /magnitude cuts. Therefore, we overestimate the total number of pairs by 62.4 
(a 2.8% correction). 

Finally, we must estimate the total number of isolated LRGs (i.e., those with no neigh- 
boring LRGs passing our CiC cuts) missing from our sample due to the incompleteness of the 
spectroscopic sample. Again making use of the assumption that objects in the photometric 
sample that would fail the color/magnitude cuts of our LRG sample are uncorrelated with 
the spectroscopic sample and that the number of expected interlopers is negligible (Eqns. [3] 
and H]), we expect the total number of isolated galaxies in the photometric sample to be 

-^photo,iso -^pQ,ss(l Pgroup) ~l~ -^faili (5) 

where p group is the probability that an LRG is in a CiC group with Usat > 1- Using the 
observable Nphoto,iso, the number of objects from the photometric catalog not grouped with 
spectroscopic objects, along with the pass rate for a random set of LRG targets of 36.6%, 
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allows us to solve for the ratio Npass/Njau in the photometric sample. Given the final esti- 



mate of the group multiplicity function (see below), we find Pgroup - 
LRG contribution from the photometric sample, Npass * (1 — Pgroup) 
an increase in total objects in ou r sample of 8%, 
reported by lMasjedi et al.l (120061 ) . 



: 0.132 and the isolated 
is 2900. This produces 
in agreement with the incompleteness rate 



2.2.2. Group-finding Algorithm and Bright Star Mask Corrections 

One-halo pairs of galaxies are assigned by satisfying the criteria Ar_|_ < r±^rnax and 
Az/{1 + z) < Pmax- Galaxies are then grouped by a FoF algorithm. The challenges in this 
section are to self-consistently incorporate the photometric sample into the calculation of 
the FoF group multiplicity function, and to marginalize over tiny Bright Star Mask holes in 
the survey geometry. 



The group finding algorithm first assigns candidate pairs for each galaxy. As de- 
scribed in § 12.2.11 the redshift space criterion for a spectroscopic-photometric pair becomes 
cimin{zspec) < ciphoto < cirnax {z spec) and -23.2 < Mg^photo < "21.2. Usiug thcsc pair as- 
signments, candidate groups are formed with the FoF algorithm. Spectroscopic neighbors of 
photometric objects in the group are kept if they satisfy the Ary requirement with at least 
one spectroscopic object in the group, and photometric neighbors of photometric objects in 
the group are kept if they satisfy the redshift space criterion with at least one spectroscopic 
object in the group. We correct the group assignments for the 23 photometric objects in 
two distinct groups. We first give preference to groups in which the photometric object is 
a fiber-collision pair; otherwise, we assign the photometric object to the larger group. We 
must then recompute the membership of the group losing the photometric object. 



The last step of the group-finding algorithm accounts for Bright Star Mask and counts 
the number of isolated photometric objects, Nphoto,iso used in Eqn. [51 The Bright Star Mask 
catalog, consisting of the ra, dec, and radius for each circu lar masked region, w as generated 
from the Tycho2 catalog as described by the NYU VAGC ( iBlanton et al.ll2005l ). The Bright 
Star Mask covers 1.88% of the total survey region. Using the observed pair counts as a func- 
tion of radius, we estimate the probability per unit area of finding a 1-halo pair satisfying 
r_L < r±^max = 0.8 Mpc/h and Az/{1 + z) < Pmax as a function of r±. The result is well-fit 
with a power law: dP/dA = 0.027rJ^'°^ (Mpc//;,)"^. dP/dA is similar to the projected cor- 
relation function, except that it does not subtract the uncorrelated contribution, and it only 
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includes pairs meeting our Ary criterion. For each bright star mask within r±^rnax of an LRG, 
we evaluate the probability that an unobserved LRG resides under the Bright Star Mask by 
integrating our power law fit to dP/dA over the area of overlap between the Bright Star Mask 
and the LRG annulus satisfying r_|_ < r_|_ ,„a^.. We assign the probability of a given Bright 
Star Mask to cover an LRG as the maximum probability assigned from all the observed 
LRG neighbors of the Bright Star Mask elements. We modify the final group multiplicity 
function to account for each possible unobserved LRG, including its ability to bridge two 
previously existing groups. Therefore a single group of LRGs can contribute to many group 
multiplicity bins with weights according to the probability of a covered LRG beneath each 
Bright Star Mask intersecting a group member. This procedure adds 79.8 net nonisolated 
LRGs to the group multiplicity function; the rest of the expected 0.0188 * iVLRc, observed are 
added to the isolated LRG count. See Tabled] for the final multiplicity function along with 
the contributions of each of these corrections. 



Finally, we place an upper limit on the contribution from photometric-object only 
groups. As listed in the final column in Table [U, we find 64 photometric object pairs and a 
single group of 3 photometric objects for which the cy and Vpetro for each galaxy in the group 
fall into the 95% range of the same redshift bin of Fig. [H 30% of such pairs fall into the 
z = 0.35 — 0.36 bin, for which cy is the least discriminating, while only 7.8% of the spec- 
troscopic LRGs lie in this redshift range. Therefore, we conclude that we miss a negligible 
number of groups by ignoring photometric object only groups. 



2.2.3. Note on average ulrg 



Zehavi et al.l (120051 ) construct a model of the spectroscopic LRG sample and estimate 
Ulrg = 9.7 x lO'^ [Mpc/h)'^ for the 0.16 < z < 0.36, -23.2 < Mg < -21.2 subsample. 
This is in good agreement with our estimate Nljig/V = 1.0 x 10~^ (Mpc//i)^^ where V 
is the product of the fraction of the sky covered by our sample and the comoving volume 
between 2; = 0.16 and z = 0.36. As discussed in the previous sections, the number of LRGs 
in our sample is increased by 8% when including objects from the imaging sample and cov- 
ered by Bright Star Masks. For our choice of cosmological para maters, the volume of the 
0.16 < z < 0.36 shell is 5% larger than in the IZehavi et al.l (120051 ) cosmology, so a difference 
of 3% in Ulrg is expected. 
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Table 1. This table lists the CiC group multiplicity function of the SDSS LRG 
spectroscopic + photometric sample for Ar_|_ max = 0.8Mpc/h and Pmax = 0.006. Column 1 

lists the number of satellites in the group, Ugat = ngroup — 1- Column 2 contains the raw 
group multiplicity before any corrections have been made (including for double counting of 

photometric objects). Column 3 lists the final group multiplicity function after all the 
corrections of § I2.2.2[ Column 4 shows the number of photometric objects in groups with 
Usat satellite galaxies, Column 5 shows the net effect of the Bright Star Mask correction, 
and Column 6 shows the number of groups of photometric objects only with Usat satellites 
satisfying the criteria listed in § I2.2.2[ These groups are largely false detections and do not 
contribute to the counts in Column 2. The total area of the survey, excluding the boundary 
region but including the bright star mask area, is 1.549 steradians. 
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2.3. Simulations and Halo Catalogs 



In order to constrain the LRG HOD, we require mock LRG catalogs to calibrate the 
relationship between the CiC group multiplicity function and the true HOD multiplicity 
function. The mock LRG catalogs are derived from halo catalogs from a. z = 0.2 A^-body 
simulation snapshot usi ng the HOD formalism detailed in § 12.41 We use a 1024'^ particle TPM 
jBode fc Ostrikerl l2003h 1 (h-^ Gpc)^ simulation w ith particle mass M p = 6.72 x lO^^M© 
and e = 16.28 h^^ kpc described in more detail in ISehgal et al.l (120071 ). The cosmological 
parameters ar e set to (T^m, ^fe, ^ a, cr^, /i) = (0.26, 0.044, 0.74, 0.95, 0.77, 0.72). For our 
FoF catalogs (IDavis et al.lll985l ). all halos containing > 107 particles were identified using a 
linking parameter 6 = 0.2. The virial mass and radius for each halo are m easured in spheres 
with overdensity determined by the spherical top hat collapse model ( iBryan &: Norman 
19981 ). We also find the radius r^^^^ at which th e circular velocity is m aximum in order 
to estimate the halo concentration Cmr- Eqn. 11 in lLokas &: MamonI (120011 ) implies 



2. 16rj;j^. /r^ 



(6) 



We produ ce catalogs u s ing t he spherical overdensity (SO) halo finder with A = 200pfe de- 



scribed in Tinker et a 



(120081) and code kindly provided by J. Tinker. We record halos down 



to 50 particles. iReidI (120081 ) shows that the clustering and velocity statistics of halos contain- 



ing 50-80 particles are not altered by the mass resolution. The mass function in this mass 
range is overestimated by ~ 5% and unaffected in higher mass bins. As we demonstrate in 
later sections, only a small fraction of LRGs occupy halos in the lowest mass bin for our best 
fit HOD, so the effect of simulation resolution on the HOD parameters should be minimal. 



2.4. HOD model 



The Halo Occupation Distribution (HOD) model assumes that the probability P^NlugIM) 
of Nlrg LRGs occupying a dark matter halo of m ass M at redshift z depends only on the 
halo mass (for a review, see ICooray fc ShethI (120021 )). However, the application of the HOD 
formalism requires that we make several further assumptions about P{Nlrg\M). Detailed 
studies of dark r natter halos and subh alos suggest a division of galaxies into central and 
satellite galaxies (IKravtsov et al.ll2004l ). The central galaxies are assumed to sit at the halo 
center, consistent with the observation that most (~ 80%) of the brightest cluster LRGs are 
found within 0.2r^,jr of the center of the cluster potential well as traced by X-rays. Satellite 
galaxies occur in the more massive halos already containing a central galax y. In high resolu- 
tion simulations they can be directly associated with dark matter subhalos (jVale fc Ostriker 
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20061 ): here we will assume they have the same distribution as the halo dark matter. We will 
use these functional forms with the five free parameters Mmm, o'logM, Mi, M^ut, and a to 
describe the mass dependence of the average halo occupation as a function of halo mass M: 



{Ncen) 



{N{M)) = (iV„ 

logioM- 



1 + erf 



„) (1 + {N,at)) 
loQioMrm 



(Nsat) 



O'logM 

' M-M,ut 

Ml 



(7) 



(9) 



T his form for N^cniM ) was suggested in IZheng et al. and a dopte d in the analysis 

of iBlake et aP J2008h . Both Isiake et all J^OOsI) and iKulkarni et all J2007h use Eqn. Iwith 
Mcut = 0. As shown below, we find evidence for M^ut > 0. Half of halos with mass Mmin 
host central LRGs, and aiogM quantifies the width of the transition from Ncen{M) = to 
Ncen{M) = 1. Ml sets the mass scale at which satellite galaxies become probable, and 
Mcut sets a cut-off below which halos do not host satellites. In the case M ^ M^ut and 
Ml ^ Mmini we expect a = 1, or the number of satellite galaxies to be proportional to the 
halo mass. 



P{NLRG,sat\N{M)) is assumed to be Poisson distr ibuted; this assumpti on is supported 
by the dist ribution of subhalo counts in s imulations (IKravtsov et al.l |2004| ) as well as ob- 
servations (ILin et al.l l2004l : IHo et al.l 120071 ). Often two-point statistics are used to fit the 
HOD; they depend only on the first and second moments of P{Nlrg\M). The technique we 
present here to constrain the HOD, Counts-In-Cylinders, makes use of higher-order statistics 
in the galaxy distribution. This technique has excellent constraining power for the param- 
eters of Nsat, but is dependent upon the accuracy of the Poisson assumption for deriving 
accurate HOD parameters. Using the Poisson assumption, the expected number of halos 
with Nsat = n satellites is given by 



{N{Nsat = n))= / dM rihaUM) exp{-Nsat{M)N,en{M)) 



{Nsat{M)N,UM)y 



(10) 



Eqn. [in] is central to our maximum likelihood analysis described in § 12.61 



The CiC method constrains only the HOD parameters in Nsat{M). In this work we 
fix the Ncen parameters M^m and aiogM by matching the observed n and amplitude of the 
projected correlation function in the 2-halo regime. 
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2.4- 1- Populating the Simulations 

Halos are populated with a central galaxy with probability Ncen{M). Central galaxies 
are placed at the center of their host halos and assigned the peculiar velocity of their halos. 
Halos with a central galaxy are populated with Ngat galaxies, where P{Nsat\N{M)) is drawn 
from a Poisson distribution. Our parameter constraints are derived using SO halo catalogs. 
For those, the position and velocity of the satellite galaxies are taken to be that of a randomly 
selected dark matter particle halo member. For mock catalogs based on the FoF halos, the 
satellite galaxies are independently distributed following an NFW profile with concentration 
of the dark matter halo determined by Eqn. El The peculiar velocity of a satellite galaxy is the 
sum of the halo peculiar velocity and a rando m velocity drawn from a Gaussian distribution 



determined by the virial velocity of the halo (ILokas &: Mamonll200ll ) 



2 _ CrMjjir (^^\ 

'^vir,lD — "^"B ■ \^^) 

We assign comoving redshift space position s to an object in our mock catalogs using the 
conversion at Zbox = 0.2: 

S = XloS + (1 + Zkox)Vp/ H{Zbox) (12) 

where xlos is the comoving distance along the line of sight in real space. 



2.5. Counts-In-Cylinders Technique 

The hypothesis underlying the Counts-In-Cylinders technique to constrain the LRG 
HOD is that 1-halo and 2-halo LRG pairs are separable based on their relative angular and 
redshift space positions. In the regime of small separations where the 1-halo term dominates, 
a cylinder should be a good approximation to the density contours surrounding central galax- 
ies, as long as the satellite velocity is uncorrelated with its distance from the halo center, and 
the relative velocity dominates the separation of central and satellite objects in the redshift 
direction. Based on our initial analysis of completeness and contamination of mock catalogs 
derived from FoF halos, we set Ar±^rnax = 0.8 Mpc/h and Azmax = 20 Mpc/h for our z = 0.2 
catalogs. Ar±^rnax is set by the typical comoving size of halos hosting satellite galaxies, and 
Azmax is set by the amplitude of the velocity dispersion in halos massive enough to host 
satellite galaxies. In later work we plan to improve the fidelity of our CiC group identifica- 
tion. However, the choice made here is sufficient since we calibrate the relatively small bias 
in the method using mock catalogs. 
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Finally, we must select a redshift dependence for the parameters Ar±^rnax and Ar^^^max- 
The virial radius of a halo of fixed mass in comoving coordinates decreases by < 6% over 
our SDSS sample redshift range; this decrease will be offset as massive halos hosting LRGs 
accrete and grow. We approximate the net result by fixing the transverse separation in 
comoving coordinates. The virial velocity of a halo of fixed mass varies by a similar fraction, 
but again we expect it to grow with halo mass. We keep the effective maximum relative 
velocity fixed as we translate the cylinder parameters to different redshifts. For an object at 
Xo = xi^d), the observed redshift will be 

Zobs = Zd+ {1 + ZdjVp/ c (13) 

so for pairs in massive halos where Azd is relatively small, Azobs/{^ + Zobs) ~ Avp/c. 



2.6. Maximum Likelihood Parameter Estimation 

While the main purpose of this technique is to produce mock catalogs with higher-order 
statistics in agreement with the observed sample (specifically the CiC multiplicity function), 
in the process we derive constraints on HOD parameters. These results are dependent on the 
input mass function, which we take directly from our simulation. The HOD model we are 
constraining predicts the expected number of groups containing Nsat satellite LRGs for each 
positive integer n = Ngat via Eqn. [TOl The CiC technique presented here produces a group 
multiplicity function which may have both an offset and scatter around the true number of 
groups with Ugat satellites, 

Ag{nsat) = Ncic{nsat) - Ntrue{nsat), (14) 

where by 'true' groups we mean groups of LRGs occupying the same halo. The overall offset 
reflects both the incompleteness of CiC method for finding pairs of LRGs in the same halo, 
as well as the contamination from LRGs in other halos. The scatter Ag{nsat) arises both 
from the stochasticity of the Nl^g ~ ^haio relation and cosmic variance of halo alignments 
that cause pairs of galaxies in separate halos to be grouped together by CiC. Using our mock 
catalogs we have estimated P{Ag{nsat)) for each Usat] see § 13.21 for more details. 

While in principle Ag{nsat) should be measured at each HOD point, this is not feasible 
since one must make many mock catalogs in order to estimate the distribution of Ag{nsat)- 
As we show below, the tightest parameter constraints from our CiC group multiplicity mea- 
surement are on Mi and a + Mc„j/(1O^^M0). Examination of Agirisat) over a subset of 
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points around the maximum likelihood HOD indicates that the variance of Ag{nsat) is much 
larger than any systematic variation of Ag{nsat) with HOD parameters for Ugat = 2, 3. There 
may be slight trends at larger Ugat- This may be because the CiC parameters adopted in 
this paper miss a significant fraction of large groups, so that a better model for the offset 
between Mac and Ntme would be Nacinsat) = Ntmeinsat) * h{nsat)- We chose not to adopt 
this parametrization because at large nsat-, either multiplicity function can be and lead to 
ambiguities in hijisat)- We hope to improve the accuracy of the CiC technique at ngat > 2 
in future work. However, we expect this effect to have only a small impact on the final 
parameter constraints. 



We assume that for each value of Usat, the number of groups with Ugat satellites, j = 
Ncicinsat) — A.g{nsat), is Poisson-distributed about the expectation value fj,n in Eqn. [101 

P{j = Nacinsat) - Agirisatm = e-'^"4 (1^) 

r- 

and independent of the other observed and expected multiplicities. Next we marginalize over 
the distribution of Ag{nsat) which we measured from our mock catalogs: 

P{Nac{nsat)W) = J P{j = Nacinsat) - Aginsat)W)PiAg)dAg. (16) 

Then the probability of observing the set of CiC multiplicities Ncic = {Ni, N2, ...Nk^^^} 
given the expectation values /2 = {fii, ^2, ••■/^fc™^} (with the subscript denoting the number 
of satellites in the group) is 

PiNcHclfl) = Ut.TPiNnlfln) (17) 

In order to make use of Bayes' Theorem to find the parameters p that maximize the likelihood 
of the observations NciC: 

Pip\Nc^c) = PiNc^c\p)Pip)/PiNc^c) (18) 

we must assume a prior Pip). We choose a flat prior on a, Mi, and Mcut, except that they 
are required to be positive. 

Finally, we note that the NacsDssinsat) = for Usat > 9. In Eqn. [T7]we set k^ax = 18. 
The exact cutoff is unimportant since nsat < 9 multiplicities already constrain the HOD 
parameter space to one in which Nirigat) is sharply falling. 
In Figures [5] through [7] we compute the expectation of several quantities as 

(y) = Jyip)Pip\Nac)dp (19) 
Ax = ^{(x-{x))^). (20) 
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Eqn. [20] is used to compute the error bars in those figures. Confidence intervals for some func- 
tion y{p) of the HOD parameters given in § 13.31 are computed directly from the marginalized 
distribution: 

P{y')dy'= [ P{p\Nc^c)dP- (21) 

Jyip)ely',y'+dy'] 

3. Results 

We use the SDSS LRG Ncic to evaluate the likelihood of a given HOD model using 
Eqn. dHl However, Mmm and <JiogM are still unconstrained. We choose (JiogM = 0.7 for the 
analysis that follows, and set so that the number density of LRGs in the mock catalog 

matches that of our SDSS sample, uljig = 1-0 x 10"'^ (Mpc/h)^^. This choice of aiogM = 0.7 
provides excellent agreement with the observed projected correlation function; see § 13.41 and 
Figure [3 

3.1. SDSS LRG Ncic 

Column 3 of Table [1] presents our final estimation of NciciP'sat) for Usat = — 8, the 
group multiplicity function for our 0.16 < z < 0.36, —23.2 < Mg < —21.2 subsample of 
SDSS LRGs. 

3.2. CiC offset Ag{nsat) and Group-Finding Accuracy 

To characterize the offsets Ag{nsat) (Eqn. UM between the CiC group multiplicity func- 
tion and the true multiplicity function, we produce 600 mock catalogs using the HOD pa- 
rameters aiogM = 0.7, M^in = 8.05 x lO^^M©, M^^t = 4.66 x IO^^Mq, Mi = 4.95 x lO^^M©, 
and a = 1.07, which are close to the final maximum likelihood parameters reported in § 13. 3j 
the two produce very similar CiC group multiplicity functions. We measured Ag{nsat) in a 
randomly selected cubic subsample of each mock catalog, selected to have a volume equal to 
that of our SDSS LRG subsample, 0.46 {Gpc/h)^. The width of P{Ag) may be underesti- 
mated, since each subsample is drawn from the same dark matter simulation. However, the 
broadness of Ncen{M) means that our simulation box contains many more halos than LRGs, 
and so the variance resulting from cosmic structures is somewhat sampled. The average 
and variance of Ag{nsat) are reported in Table El Comparison to the observed Ncicinsat) 
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in Table [T] shows that the width of the Ag distribution produces comparable uncertainty in 
{Ntrue{nsat)) as the Poisson sampling. Therefore, neglecting the integral in Eqn. [T6] would 
cause an underestimation of the errors on the HOD parameters. In our likelihood calculation 
we use the histogram of Ag^Ugat) values to estimate P{Ag{nsat)) in Eqn. [T^ 

Table [2] also reports a measure of the completeness and contamination with which we 
find groups of size Usat- Column 5 is the fraction of true one-halo groups with ngat satellites 
that are exact matches to a CiC group. This is a rather stringent definition, since it excludes 
groups with only one missing satellite of many, or with only one galaxy from a different 
halo. Column 6 shows the fraction of CiC groups that do not exactly match a true one-halo 
group. Our optimization mainly focused on groups with only one satellite, since they are 
most numerous. Moreover, the CiC parameters were established while still using FoF halo 
catalogs; there are non-negligible differences in the two-halo term between the two catalogs 
(see §3.5p . Therefore, we are still optimistic that we can improve the accuracy of our method 
at Usat > 1, which we will address in a later paper. The issue is of less concern in the current 
work, since we have calibrated the observed statistics on our mock catalogs. Large biases are 
unlikely as long as the mock catalogs are reproducing the salient features of the small-scale 
clustering and distribution of satellites within dark matter halos. 



3.3. HOD Constraints 

With the width of the Ncen{M) distribution fixed at aiogM = 0.7, we compute the likeli- 
hood at each point in the three dimensional space of N^at parameters a, Mcut, and Mi accord- 
ing to Eqns. [15]- [181 We evaluate Mmin at each HOD point so that ulfig remains fixed. The 
maximum likelihood HOD parameters and marginalized one- dimensional 68% and 95% con- 
fidence intervals are M^^t = 5.0lJJ(t|^) x W^Mq, Mi = 4:.95toil{toll) x IO^^Mq, and a = 
IMbtoijitoit)^ with M^in = 8.05 x W^Mq at the maximum likelihood point. In Figures[21 
|3l and IHwe show two dimensional likelihood contours for AL = {—1.15, —3.09, —5.9} after 
marginalizing over the remaining parameter. For a distribution, these contour values cor- 
respond to 1,2, and 3a confidence regions. We find a strong degeneracy in the Mcut — en plane. 
We tightly constrain the fraction of galaxies that are satellites, fsat = 0.0636lQ;oo2o(-aoo39) 
and the parameter combination a + Mcut / i^O^'^ Mq) = ^■^'^-oMri-omo)- orthogonal 
direction, we find only a weak constraint: Mcut/ {10^'^ Mq) - a = -O.bAtoHitoAg)- In Col- 
umn 3 of Table [3] we list the number of groups with Ugat satellites averaged over 20 mock 
catalogs evaluated at our maximum likelihood HOD, which show good agreement with the 
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Table 2. Column 2 shows the average offset Ag^Ugat) = Ncicingat) — Ntrueingat) over 600 
mock catalogs for HOD parameters aiogu = 0.7, M^m = 8.05 x lO^^M©, 
Mcut = 4.66 X IO^^Mq, Ml = 4.95 x 10^^ Mq, and a = 1.07, which are close to the final 
maximum likelihood parameters reported in § 13.31 Column 3, a^, is our estimate of the 
variance of Ag in a subsample of the simulation box with volume equal to our SDSS 
subsample, 0.46 [Gpc/h)^. Comparison with Table [1], Column 2 (reproduced here in 
Column 4) shows that this variance is comparable to the expected variance from Poisson 
sampling of {Ncicingat)) ■ Column 5 is the fraction of true one-halo groups with Ugat 
satellites that are exact matches to a CiC group. Column 6 shows the fraction of CiC 
groups that do not exactly match a true one-halo group. 
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SDSS CiC multplicity function; recall that the broad distribution in Ncic — Ngat introduces 
noise in addition to the Poisson term. 



In Figures [5] through [TJ we show the implications of these results for the distribution 
of both satellite galaxies as a whole and groups containing Usat satellites as a function of 
halo mass. Figure [5] shows the mean and rms (Eqns. [19] and [20l) of the expected number 
of satellites, {Ncen{M)Nsat{M)) for several halo masses. Nsat{M) is less well constrained at 
the high mass end, where the number density of halos is low. Furthermore, we have not 
included cosmic variance of halo counts in this analysis; we have only used the mass function 
from our simulation. By comparison with Figure [3, we see that Ncen{M) < 1 even when M 
is large enough to host satellite galaxies. Thus the tight constraints on the low mass end of 
Nsat{M) are dependent on the accuracy of our parametrization of NceniM). Figure [6] shows 
the probability that a halo of mass M hosts n^at satellites. Again, the constraints are weaker 
as M increases and n{M) decreases. Figure [7] shows the distribution of all satellite galaxies 
as a function of halo mass, as well as the distribution of satellites in groups with n^at = 1,2,3 
and 4. The width of these distributions is comparable to the difference in mean halo mass 
as function of n,^at ; this is qualitatively in line with the large scatter seen in N^fiG vs. M200 



m 



Ho et al.l (I2OO7I ). Attempts to measure the mass of dark matter halos using LRG groups 



as a tracer should expect a broad distribution. 
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Fig. 2.— Contours for AlnL = {-1.15,-3.09,-5.9} in the Mi vs. M^ut plane after 
marginalizing over a. For a distribution, these contours would enclose 1, 2, and 3a" 
confidence regions. The cross indicates the maximum likelihood parameter values. 
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Fig. 5. — The dashed curve ranging from to 1 shows the Ncen{M) term for the maximum 
hkehhood HOD; it should vary only slightly with HOD since the satellite fraction is well 
constrained by our model, and we hold criogM fixed. The solid curve shows {Nsat{M)Ncen{M)), 
with the error bars computed by Eqn. [201 
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Fig. 6. — The probability Pijisat = MM) (long-short dashed curve), Pijisat = 2|M) (long 
dashed curve), Pijisat = 3|M) (short dashed curve), P^rigat = 4|M) (solid curve) vs. halo 
mass M/Msun- Error bars are computed by Eqn. [201 Though P{nsat\M) are evaluated at 
identical values of M, points are slightly staggered for clarity. 
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Fig. 7. — The total number density of satellite galaxies occupying a halo of mass M per 
IoqiqM (thick solid upper curve). The long-short dashed curve shows the number density 
of satellites occupying halos with Ugat = 1, long dashed curve for satellites in halos with 
Usat = 2, short dashed curve for ngat = 3, and solid curve for Ugat = 4. Error bars are 
computed by Eqn. 



- 28 - 



3.4. Projected Correlation Function Wp 



Zehavi et al.l (120051 ) and lMasjedi et al.l (120061 ) measure the projected correlation function 



Wp{rp) for this sample on small and intermediate scales: 



Wp{rp) = 2 rf<(rp,7r) 
Jo 



(22) 



We follow IZehavi et al.l (120051 ) and set vr, 



80 Mpc/h, whi ch is large enough to include 



most correlated pairs and produce stable estimates of Wp{rp). iMasjedi et al.l (120061) recover 
missing fiber collision pairs by computing Wp{rp) by cross correlation between the SDSS 
spectroscopic and imaging samples. They also correct for photometric biases of close galaxy 
pairs, which can introduce incompleteness of pairs with separation < 0.1 Mpc//i. 



We present the projected correlation function Wp{rp) averaged over 20 mock catalogs 
produced with our SO halo catalog using our ma ximum likelihood HO D in Figure [H We 
find excellent agreeme nt with the measurem ents of lMasjedi et al.l (120061 ). Using the diagonal 
error bars reported in lMasjedi et al.l (120061 ). we find = 7.5 for the outer 15 points. There 
is substantial discrepancy with the inner 3 points at Vp = 0.01,0.016,0.026 (not shown in 
Fig. E]); = 29 for all 18 points. The discrepancy is not surprising since these small 
distances are comparable to the force resolution of our simulation. Though the CiC method 
relies primarily on pairs with Vp < 0.8 Mpc/h, our mock catalogs reproduce the features of 
the observed Wp{rp) by adjusting a single parameter aiogM to match the large scale 



Mpc//i) bias probed by w. 



20 

Note that a sharp transition from to 1 for Ncen{M) can be 
ruled out with confidence. Figure [8] shows Wp{rp) for catalogs with aiogM = 0.2, 0.7, and 1.3 
for comparison. All three catalogs match the observed clustering at Vp < 0.8 Mpc/h where 
we have CiC constraints, but only catalogs with aiogM ~ 0.7 match the observed clustering 
on ~ 2 — 20 Mpc/h scales, the regime where two-halo pairs dominate. 



- 29 - 



1000 



100 



10 




Fig. 8. — The projected correlation function Wp(rp) vs. pro j ected separation r^, both in 



Mpc/h. The points with error bars taken from lMasjedi et aLl (120061 ). The sohd curve is an 
average over 20 mock cat alogs at the maximum likelihood HOD parameters presented in the 
text for aiogM = 0.7. The lMasjedi et al.l (120061 ) data was used in our analysis only to optimize 
the value of criogM, which sets the large scale clustering. The striking agreement shows that 
using SO halo catalogs we can simultaneously reproduce the CiC group multiplicity function 
and the clustering on scales ~ 1 — 20 Mpc/h. For comparison, we produce a catalog with 
<^iogM = 0.2 (long dashed curve) and ciogM = 1-3 (short dashed curve). Both catalogs match 
the observed clustering at Vp < 0.8 Mpc/h where we have CiC constraints, but disagree with 
the observations on larger scales. 
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3.5. Comparison with Results using FoF Catalogs and Previous Works 



We also generate mock LRG catalogs based on a FoF halo cat alog produced with linking 
length b = 0.2. For this catalog we use virial masses as defined in lBryan Sz Norman 
Setting aiogM = 0.5, we are able to match both the SDSS multiplicity function (Column 4 
of Table [3]) and the large scale projected correlation. However, Figure [3 shows that the FoF 
catalog does not repro duce the observed w Jrp) at rp ~ 1 Mpc//i. Using the diagonal error 
bars only presented in lMasjedi et al.l (120061 ). this catalog has = 104 for the last 15 points 
(the mock catalog based on SO halos has = 7.5). The large discrepancy between the 
FoF and SO-based mock catalogs arise s from the FoF algo rithm's tendency to link nearby 
halos. The SO algorithm explicated in [Tinker et al.l (120081 ) allows halo overlap in order to 
correctly separate distinct density peaks likely to host LRGs. This problem should also be 
circumvented by algorithms which identify subhalos. 



If one sums the difference between the observed and FoF catalog Wp{rp) for the three 
points near 1 Mpc/h, the number of missing pairs in that regime is equal to half the number 
in the one-halo pair regime, Vp < 0.8 Mpc/h. To match only Wp{rp) using an FoF catalog, 
one would therefore need to enhance the number of pairs at separation ~ 1 Mpc/h without 
greatly altering the number of pairs at smaller radii. In the halo model, this can be accom- 
plished by putting many satellite galaxies in a large halos whose virial radii extend to ~ 1 
Mpc/h, so that satellite-satellite pairs also contribute substantially to Wp{rp). In a given halo, 
satellite-satellite pairs will have larger separations than central-satellite pairs. This drives 
a to large values so that the one-halo term can accommodate the missing two-halo pairs in 
FoF catalogs. Current LRG HOD results in the literature may theref ore be systematically 
biased if the two-halo term is derived from F oF halo catalogs directly (IKulkarni et al.l 120071 : 
White et al.l 120071 : iPadmanabhan et al.ll2008l). based on ana l ytic fits to the HOD calibrated 



on Fo F catalogs from A^-body simulations (iBlake et al.ll2008l: iBrown et al.ll2008l : IZheng et al. 
20081 ). or neglects the scale-dependent bias (IWake et al.ll2008l ). 



Both IZheng et al.l (120081 ) and iKulkarni et al.l (120071 ) report HOD constraints for the 
SDSS LRG subsample studied in this paper. Our tight constraint on the satellite f r action , 
6.4 ± 0.4%, is in excellent agreement with the 5.2 — 6.2% found by IZheng et al.l (120081 ). 
However, they find a steeper value of a, 1.8 at erg = 0-8) than our finding of a = l-Olg g 
(95% confidence). Equivalently, they find {NceniM)Nsat{M)) ~ 4 at M = IO^^Mq while we 
find ~ 2. This difference is p robably due to the diff erence in our assumed 2-halo correlation 
function at ~ Mpc/h scales. IKulkarni et al.l (120071 ) provide their halo mass function, so we 
are able to make a direct comparison between their predicted multiplicity function and our 
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Ncic measurement. Column 5 of Table [3] shows the expected true multiplicity function from 
the best fit HOD presented in lKulkarni et al.l (120071 ). scaled to an SDSS volume. Their HOD 
predicts many more groups with Ugat = 1,2, and 3 than are observed; at larger n^at our 
current CiC method may not recover such large groups. This HOD also predicts 15 objects 
in an SDSS volume with Ugat > 15. Interestingly, if we fix Mcut = 0, we also find a ~ 1.4 
(see Figure [3]). The cause of the discrepancy here seems to be mainly in Mi, the overall 
normalization of their Nsat{M). While we find fsat = 0.064, they report fsat = 0.17. Per- 
haps their large FoF linking length b = 0.6 Mpc/h severely reduces the two-halo clustering 
in small scales, forcing the one-halo term to compensate. 



The results of iHo et al.l (120071 ) su pport our findin gs for a relatively low cutoff for the 



number of LRGs per halo. Figure 12 of iHo et al.l (120071 ) shows a single cluster with 16 LRGs; 
and all other clusters in their sample have < 11 LRGs. We expect this to provide a strict 
upper bound on our sample, since their photometric sample has a number density ~ 4 times 
larger than the spectroscopic sample we study here. 



4. Assessment of the CiC method 

The HOD parameter constraints reported in § 13.31 are for aiogM held fixed. Since the 
CiC constraints always involve the product {Ncen{M)Nsat{M)) and Figure [5] shows that 
Ncen{M) < 1 in regions where the satellite contribution is non-negligible, the CiC maximum 
likelihood HOD parameters will vary with aiogM to hold {Ncen{M)Nsat{M)) relatively fixed. 
Moreover, varying aiogM changes the average bias of the LRG-occupied halos, which may in 
turn affect the distribution of Ag{nsat) m Eqn. [HI We expect this to be a smaller effect. A 
complete HOD analysis would include Ncen{M) parameters when evaluating the likelihood. 
We find that using Wp{rp) on ~ 20 Mpc/h to set <JiogM is sufficient for the purposes of gen- 
erating mock catalogs. We hope to explore variations in the functional form of Ncen{M) 
in later work. We have not explored the dep endence of our method on the assumed distri- 
bution of satellite galaxies. ISeo et al.l (120071 ) find that passively evolving galaxies are more 



concentrated than the dark matter at low redshift; this should make the CiC group finder 
more accurate than we found for our mock catalogs. For the parameters chosen here, the 
CiC method does not extract exact matches to one-halo groups with high fidelity. We have 
circumvented this issue by careful calibration on our simulated catalogs, and we expect to 
refine the method in the future to improve the accuracy of finding groups with more than 
two LRGs. 
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Fig. 9. — The projected correlation function Wp(rp) vs. projected separ ation r^, both in 
Mpc/h. The points with error bars from taken from iMasjedi et aLl (120061 ). The sohd curve 
is an average over 20 mock catalogs at the maximum likelihood HOD parameters presented 
in the text for aiogM = 0.7 and using the SO halo catalog. The short dashed curve shows 
Wp{rp) averaged over 5 mock catalogs generated by populating an FoF halo catalog. The 
FoF halo catalog has a severe suppression of two-halo pairs at separations of ~ 1 Mpc/h, 
equal to the half of the total number of pairs in the one-halo regime, Vp < 0.8 Mpc/h. 
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Because the CiC method can simultaneously reproduce the large scale clustering of 
LRGs and match the higher order statistics in the density field probed by the CiC group 
multiplicity function, we expect it to be a useful tool both for HOD constraints and the 
production of good mock catalogs. Since the CiC method only uses information on scales 
where the one-halo term dominates, our method is less sensitive to knowing the two-halo 
correlation function on ~ Mpc//i scales. We have demonstrated that the deficit of two-halo 
pairs at ~ 1 Mpc/h is problematic for HOD fitting methods using 2 and 3-pt statistics and 
FoF catalogs because the HOD must accommodate the missing two halo pairs at this sepa- 
ration, and thus predict groups containing many more LRGs than are observed. 

Finally, the low number density of the LRGs makes them particularly well-suited for 
this method. 



5. Conclusions 

The low number density of SDSS LRGs allows us to partially separate LRG pairs occu- 
pying the same dark matter halo from pairs occupying distinct dark matter halos. Candidate 
one-halo pairs are identified using simple cuts in the transverse and LOS separations. We 
group these pairs using the FoF algorithm to compute the CiC group multiplicity function. 
We measure the CiC group multiphcity function for the subsample of SDSS LRGs satisfying 
—23.2 < Mg < —21.2 and 0.16 < z < 0.36, carefully accounting for the effects of fiber 
collisions and survey boundaries, holes, and incompleteness. 

In order to derive HOD constraints from our measurement, we calibrated the relation 
between the CiC and true one-halo group multiplicity functions using mock LRG catalogs. 
The variance about the mean relation is comparable to the Poisson sampling variance of the 
CiC multiplicity function and must be properly accounted for in the maximum likelihood 
parameter estimation. 

The CiC group multiplicity function places strong constraints on the satellite LRG HOD, 
NsatiM). When we fix aiogM — 0.7 and Ulrg — 10~^ (Mpc//i)~^, the maximum likelihood 
HOD parameters and marginalized one-dimensional 68% and 95% confidence intervals are 
M,„, = b.Otliitli) X lO^^Me, M, = 4.95lg|,^(lLl) x IO^^Mq, and a = 1.035^gi?(lgit), 
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with Mr, 



8.05 X IO^^Mq at the maximum likelihood point. We tightly constrain the 



satellite fraction to fsat = 0.0636_Q;QQ2o(-aoo39)- '^^^ projected correlation function Wp{rp) 
of mock cata logs derived fr o m an SO h alo catalog is an ex cellent agreement with the mea- 
surements of iMasjedi et al.l (120061 ) and IZehavi et al.l (120051 ) when the large scale clustering 
is used to fix aiogM- Fig. [9] shows that FoF halo catalogs have a severe deficit of pairs at 
~ 1 Mpc//i. In § 13.51 we point out that methods using Wpivp) with an analytic estimate of 
the two-halo term calibrated using FoF halos will severely overestimate the number of satel- 
lites and the maximum expected one-halo group size. Our mea sured CiC group multiplicity 
function rules out the best fit HOD from iKulkarni et al.l (120071 ). 



Despite the increased complexity of our approach and necessary calibration using sim- 
ulations, we have produced high quality mock catalogs which reproduce both higher order 
statistics in the density field and the features of the projected correlation function. These 
mock catalogs will be used in a forthcoming paper to study the large scale structure statistics 
of our CiC groups. 
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Table 3. Column 1 is the number of satellites in the group. Column 2 is a copy of Column 
3 of Table [H our best estimate of the SDSS CiC multiphcity function. Column 3 shows the 
CiC multiplicity function averaged over 20 mock catalogs produced with our maximum 
likelihood HOD using the SO halo catalogs and scaled by the volume ratio VsDssiysim- 
Column 4 shows the CiC multiplicity function for mocks made wi th the FoF cata l og. 
Column 5 shows our estimate of the expected Ntme in Vsdss for the iKulkarni et al.l (120071 ) 

HOD. 
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